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I THE  PROBLEM 

In  (2)  we  have  discussed  the  importance  of  techniques  for  predicting 
success  differentially  in  each  of  a number  of  different  activities  from  a 
single  battery  of  predictors.  It  was  assumed  that  intercorrelations  for  a 
large  battery  of  predictor  variables  were  available  and  also  correlations 
between  these  predictors  and  a large  nxmiber  of  criterion  variables.  The 
pi*oblem  was  to  select  from  this  larger  battery  of  predictors  that  subset  of 
specified  size  which  would  yield  the  maxijnum  index  of  differential  prediction 
for  the  criterion  variables.  The  index  of  differential  prediction  efficiency 
was  taken  to  be  a siiiiple  function  of  the  average  of  the  variances  for  the 
predicted  difference  scores  for  all  possible  pairs  of  criterion  variables. 

The  larger  this  average  variance  the  greater  the  differential  prediction 
efficiency  of  the  battery.  It  was  sho\m  that  this  index  is  equivalent  to  the 
difference  between  the  average  variance  of  the  predicted  criterion  measu/'es 
and  the  average  of  their  covariances,  assuming  standard  measures  for  both 
predictors  and  criteria  and  that  the  predicted  criteria  are  the  "least  square" 
estimates.  A method  for  selecting  that  subset  of  predictors  of  specified 
size  which  would  yield  the  maximum  index  of  differential  prediction  was 
presented. 

The  method  referred  to  tacitly  assumes  that  all  predictors  in  the  battery 
take  the  same  amount  of  administration  time,  so  that  all  subsets  of  the  sams 
size  would  also  take  the  same  amount  of  administration  time.  Usually  this 
will  not  be  the  case.  A more  general  AT"iroach  to  the  problem  might  be  to 
start  with  a given  battery  of  predictor  variables  and  inquire  how  the  admin- 
istration time  for  each  predictor  should  be  altered  so  that  for  a specified 
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over-all  testing  time  the  index  of  differential  prediction  efficiency  will 
be  a maximum.  This  approach  wo’iLd  allow  for  increasing  the  length  of  an 
experimental  battery  as  well  as  for  decreasing  it. 

As  a matter  of  fact  fir  the  case  of  a single  criterion  a method  is 
already  available  (1)  for  determining  the  optimal  distribution  of  testing 
time  for  a battery  of  predictors,  assuming  that  intercorrelation,  validity 
and  reliability  data  are  available  for  predictors  of  arbitrary  lengths.  It 
is  the  purpose  of  this  article  to  present  a modification  and  generalization 
of  the  method  for  the  case  of  differential  prediction  involving  a ntunber  of 
criterion  variables. 

In  this  presentation  testing  time  is  taken  to  be  the  time  actually 
allotted  the  examinee  for  taking  the  test.  A more  complete  analysis  must 
also  take  into  account  the  time  for  reading  instructions,  practice  exercises, 
passing  out  and  collecting  papers,  etc.  The  method  will  first  be  described 
and  illustrated  by  a numerical  example  after  which  the  mathematical  rationale 
will  be  presented. 


II  HUMERICAL  EXA^iPI£ 

The  predictor  variables  used  in  this  example  are: 

(1)  Guilford-Zinnierman  Aptitude  Survey,  Part  I,  Verbal  Coi^^rehension 

(2)  Ckiili  ^rd-Juiuorman  .ptltudo  SujrVvy,  Part  HI,  'umerical  Operations 

(3)  Guilford-Ziianerman  Aptitude  Survey,  Part  VII,  Mechanical  Kno^viedge 

(4)  A.  C.  E.  Psychological  Examination,  Quantitative  Reasoning 

(5)  A.  C.  E.  Psychological  Examination,  Linguistic  Reasoning 

(6)  Cooperative  English  Test  (Form  OM),  Usage 

The  matrix  of  test  intercorrelations  with  reliabili cies  in  the  diagonal 
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is  given  in  Table  1, 


Table  1 


R Matrix  of  Predictor  Intercorrelations  v/ith  Reliabilities 
Substituted  for  Unities  in  the  Diagonal: 


R 


r-D 


1 2 3 '‘  I " T'_  6 g- 


1 

G-Z  1 

.920 

.159 

.152 

.281 

.763 

.515 

2.790 

2 

G-Z  3 

.159 

.920 

.003 

.369 

.292 

.243 

1.986 

G-Z  7 

.152 

.003 

.920 

.200 

.142 

-.150 

1.267 

4 

ACE-Q 

.281 

.369 

.200 

.820 

.549 

,426 

2.645 

5 

ACE-L 

.763 

.292 

,142 

.549 

.330 

.628 

3.2O/+ 

6 

English 

.515 

.243 

-.150 

.426 

.628 

.860 

2.522 

L 

2.790 

1.986 

1.267 

2,645 

3.204 

2.522 

14.414 

The  criterion  variables  are  grade  point  averages  in  each  of  ten  college 
subjects.  The  matrix  of  validity  coefficients  is  given  in  Table  2, 


Table  2 

The  r^  Matrix  of  Validity  Coefficients 


1 

G-Z  1 

2 

G-Z  3 

"3"" 

G-Z  7 

4 

ACE-Q. 

5 

ACE-L 

English 

E 

1 Anthropology 

.370 

.177 

.091 

.294 

.341 

.357 

1.630 

2 Chemistry 

.317 

.274 

.016 

.309 

.364 

.399 

1.679 

3 Economics 

.339 

.2U 

.008 

.241 

.334 

.323 

1.456 

4 English 

.526 

.247 

-.075 

.262 

.488 

.524 

1.972 

5 Foreign  Lang. 

,295 

.267 

-.156 

.200 

.232 

.426 

1.284 

6 Geology 

.184 

.140 

.094 

.170 

.229 

.214 

1.031 

7 History 

.379 

.169 

001 

.182 

.373 

.336 

1.438 

8 Mathematics 

.287 

.348 

—.088 

.350 

.336 

,401 

1.634 

9 Psychology 

,41-0 

.170 

.096 

.235 

.409 

.403 

I.8O3 

10  Zoology 

.336 

.a6 

.031 

.318 

.345 

.351 

1.597 

r 

3.473 

2.239 

.016 

2.611 

3.451 

3,734 

15.524 

e/10 

.347 

.224 

.002 

.261 

.345 

.373 

1,552 

The  over-all  testing  time  for  the  tests  of  arbitrary  length  is  l42 
minutes,  VJe  assume  that  this  time  is  to  be  cut  in  half  so  that  tho  over^-all 
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testing  time  is  71  minutes.  The  problem  is  to  determine  the  time  to  be 
allotted  to  each  test  so  as  to  maximize  the  index  of  differential  prediction 
efficiency 5 

The  traditional  assumptions  are  used  here  as  in  (l)  with  respect  to  the 
effect  of  test  length  on  correlation  and  vdll  not  be  repeated. 

The  method  of  solution  for  the  new  test  lengths  involves  a series  of 
successive  approximations.  For  large  numbers  of  predictor  and  criterion 
variables  the  solution  may  become  very  laborious.  It  is  probable  that  the 
solution  would  be  greatly  expedited  by  the  use  of  high  speed  computing  equip- 
ment. Further  research  may  yield  more  efficient  computational  procedures. 

1.  The  first  computational  step  is  to  calculate  a matrix  from  the 

matrix  r'  in  Table  2.  The  elements  a of  Table  3 are  the  corresponding 
c c 

elements  of  Table  2 with  their  col\amn  means  subtracted.  Hence  the  columns 
of  Table  3 all  add  to  zero. 


Table  3 

The  a'  Matrix;  Validity  Coefficients  Expressed 
in  Deviation  Form  for  Each  Test 


1 

2 

3 

4 

5 

T” 

1 

.023 

-.C47 

.089 

o' 

-.004 

-.016 

2 

-.030 

.050 

.014 

.048 

.019 

.026 

3 

-.008 

-.013 

.006 

-.020 

-.011 

-.050 

4 

.179 

.023 

-.077 

.001 

.143 

.151 

5 

-.052 

.063 

-.158 

-.061 

-.113 

.053 

6 

-.163 

-.084 

.092 

-.091 

-.116 

-.159 

7 

.032 

-.056 

-.003 

-.079 

.023 

-.037 

8 

-.060 

.124 

-.690 

.089 

-.009 

.028 

9 

.093 

-.054 

.094 

.024 

.064 

.030 

10 

-.011 

-.008 

.029 

.057 

.000 

-.022 

Ck 

.003 

-.001 

-.004 

.001 

.001 

.004 

2 

0 

0 

-.001 

-.004 

.001 

.001 

.004 

2.  T 1=  ne.'ct  step  is  to  computo  the  cleiiients  for  a diagonal  matrix,  A . 

The  i'th  element  is  the  product  of  the  original  length  of  test  i multiplied 

« 

by  one  minus  its  reliability.  The  elements  for  A are  given  in  row  4,  labeled 
I'a  of  Table  4.  For  the  first  element  we  have  A^^  » 25(1*00  - ,92)  = 2,00  . 

Table  4 

Computation  of  l^D^  pnd 


/ \ 

First  approximation:  I'Dj^  “/  )1'D 

°1  \ / 


1 

2 

3 

4 

5 

S Cic 

Z 

1 I'd 

a • 

25.0 

9.0 

30.0 

23.0 

15.0 

40.0 

142.0 

12.5 

4.5 

15.0 

11.5 

7.5 

20.0  71.0 

71.0 

3 I'DJ^ 

.080 

.222 

.0667 

.087 

.1333 

.050 

4 1'a 

2.000 

.720 

2.400 

4.140 

2.550 

5.600 

17.A10 

5 I'AD-i^ 

.160 

.160 

,160 

.360 

.340 

.280 

1.460 

3.  A first  approximation  is  now  required  for  the  altered  test  lengths.  We 

assume  the  new  test  lengths  to  be  proportional  to  the  original  test  lengths. 

Therefore,  as  a first  approximation  to  the  new  test  lengths  v;e  Lake  one  half  the 

original  test  lengths.  Row  1 in  Table  4 gives  the  crigin.'.l  test  lengths.  Row  2 

of  the  same  table  is  one  half  the  first, 

4.  Calculate  the  reciprocals  of  the  elaments.  These  are  given  in 

row  3 of  Table  4. 

5.  Calculate  the  product  of  each  A »ralue  in  row  4 of  Table  4 by  the 
corresponding  value  immediately  above  it.  The  ratios  are  entered  in  row  5 of 
Table  4.  For  example,  the  first  value  is  ,l60  = 2.000  x ,030  . 

6.  Next  the  elements  calculated  in  step  5 are  added  to  the  corresponding 
diagonal  elements  of  Table  1,  and  the  table  is  copied  into  the  up^er  left 
quadrant  of  Table  5*  The  first  diagonal  element  is  ?„,080  = .160  + .920  . 
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Note  that  the  elements  below  the  diagonal  are  not  copied  in.  The  upper  right 
section  of  Table  5 is  a , the  transpose  of  Table  3« 

w 

7.  Vfe  next  calculate  a matrix  by  preraultiplying  the  matrix  by  the 
inverse  of  the  matrix  in  the  upper  left  quadrant  of  Table  5.  The  computa- 
tions for  the  forward  solution  are  given  in  the  two  lower  quadrants  of  Table  5 
and  in  Table  6.  The  back  solution  is  given  in  Table  ?•  The  procedxire  for 
multiplying  a matrix  by  the  inverse  of  a symmetric  matrix  is  outlined  in  (3)« 

8.  The  second  approximation  to  the  new  test  lengths  is  computed  in  the 
lower  section  of  Table  7 as  follows: 

Row  a consists  of  the  sum  of  squares  of  column  elements  of  the  matrix. 
For  exainple,  the  first  element  in  row  a,  namely,  .0626,  is  the  sum  of  squares 

of  the  first  ten  elements  in  column  1 of  Table  7* 

Rovf  b is  copied  from  rov/  4 of  Table  4. 

Row  c consists  of  the  products  of  corresponding  elements  in  the  two  pre- 
ceding lines,  For  example,  .1251  - .0626  x 2.00  . 

Row  d consists  of  the  squai’e  roots  of  corresponding  entries  in  the  pre- 
ceding line.  For  example,  .3537  “ . The  computations  to  the  right 

of  this  line  and  designated  s are  obtained  by  dividing  the  over-all  new  test- 
ing time,  71  minutes,  by  1.8823,  the  sura  of  the  elements  in  the  row.  This 
gives  s = 37. 7198. 

Row  e is  a check  row.  Each  element  in  the  second  line  above  it  is 
divided  by  the  element  immediately  above.  Thus  = .3537  • 

Row  f is  obtained  by  multiplying  each  elranent  in  row  d by  s.  For  example, 

the  first  element  is  13,3415  = .3537  x 37,7198  • This  line  gives  the  second 
approximation  to  the  new  test  lengths. 

Row  g is  obtained  by  dividing  each  element  in  the  preceding  row  into  the 

2 00 

corresponding  value  in  row  b.  For  example,  the  first  value  is  .I50  = * 
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Computation  of  (R+AD^  )'^  a * L,  continued 


Computation  of  (R+AD^  )*^  a = L,  Back  Solution 
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Row  h is  a check  on  the  preceding  row.  Each  element  iii  this  ro\i  is 
multiplied  by  the  corresponding  element  in  rov/  f to  give  row  b.  For  example, 
the  firct  eleraent  is  2.001  = .I5O  x 13.3‘+15* 

Row  i is  obtained  by  adding  the  elements  in  rew  g to  the  corresponding 
reliabilities.  For  example,  the  first  element  is  i.070  = .I50  .920. 

9.  A new  matrix  is  now  computed  by  repeating  steps  6 and  7 and  using 
the  elements  of  row  i of  Table  7 in  the  diagcnal  positions  of  Table  1.  Tho 
new  matrix  is  given  in  transposed  forn  in  Table  8,  rows  1 through  10. 

10.  Step  8 is  repeated  in  rows  a through  i of  Table  8.  Row  f of  this 
table  gives  a third  approximation  to  the  altered  test  lengths. 

Steps  6,  7 and  8 are  repeated  to  get  succeeding  approximations  to  the 
test  lengths.  The  calculations  were  carried  to  5 successive  approximations  for 
the  new  test  lengths,  not  counting  the  first.  These  are  summarized  in  Table  9. 
As  will  be  seen,  the  iterations  have  not  completely  stabilized.  However,  for 
practical  purposes,  the  approximation  is  doubtless  adequate. 

11.  To  compute  the  successive  indices  of  differential  prediction 
efficiency,  0^  , we  proceed  as  follows; 

(a)  For  the  index  corresponding  to  the  first  approximation  to  the  nsw 
test  length  multiply  each  element  in  the  matrix  in  Table  7 by  the  corres- 
ponding element  of  Table  3 and  sura  the  products.  This  is  the  first  entry, 

.227,  in  the  0 column  at  the  right  of  Table  9* 

(b)  To  got  02  follow  the  same  procedure  except  use  the  matrix  in 
Table  8 instead  of  in  Table  7. 

(c)  In  the  same  way  calculate  subsequent  0's  by  using  the  elements  in 
the  corresponding  L matrix  and  the  elements  in  Table  3* 
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Table  9 

Successive  Approximations  to  l^D^ 


for 


71 


Aoprox'  n 

1 

2 

3 

— 

. 5._ 

Value 

Cl  1 A A e 

of  0 for 
3Sive  Vail 
of  L 

(.5)1'E  : 

1 

12.50 

4.50 

15.00 

11.50 

7.50 

■20.00  • 

71.00 

olXC  C 

2 

13.34 

5.50 

12.76 

12.57 

12.55 

14.28 

71.00 

.227 

3 

13.27 

5.32 

11.55 

12.52 

16.05 

12.29 

71.00 

.234 

4 

13.23 

5.20 

10.98 

12c47 

17.62 

11.51 

71 . 01 

h 

.235 

5 

13.31 

5.15 

10.76 

12.46 

18.13 

11.19 

71.00 

.236 

6 

13.35 

5.12 

10.70 

12.46 

18.37 

11.00 

71.00 

.2yi 

It  \/ill  be  noted  that  f does  not  increaso  much  in  this  particular  illus- 
tration. It  goes  from  ,227,  taking  the  test  lengths  as  one  half  their 
original  length,  to  ,237  as  they  approach  optimal  length.  This  is  an  increase 
of  less  than  5 pe^'  cent  even  though  several  of  the  test  lengths  are  changed 
greatly.  For  example,  test  5 increases  f^oin  7 to  IG  minutes  while  test  6 
reduces  from  20  to  11  minutes. 

The  technique  was  applied  to  tho  same  data  assuming  that  the  total 
administration  time  was  to  be  the  same  as  in  the  original  administration, 
namely,  142  minutes  and  also  assuming  it  was  to  be  doubled  to  284  minutes. 

Only  three  approximations  to  the  optimal  test  lengths  were  calculated  for 
each  of  these  two  conditions.  Tables  10  and  11  summarize  the  results  for 
the  two  conditions  respectively.  The  last  column  in  each  table  shows  the 
index  of  differential  prediction  efficiency,  0,  corresponding  to  each 
approximation  to  optimal  test  lengths.  In  both  cases  the  improvement  of  0 
as  tho  tests  approach  optimal  length  is  appreciably  greater  than  for  the  case 
of  one  half  the  original  testing  time.  This  rate  of  improvement  is  greatest 
for  double  testing  time.  As  can  by  seen  from  the  right  hand  column  of  Table  11 
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it  goes  fi’om  .305  to  .337  wliich  is  approxiinate]  y a 10  per  cent  increase. 
Further  research  is  needed  to  determine  the  sensitivity  of  0 to  alterations 
in  relative  testing  time  for  each  of  the  tests  and  to  variation  in  total 
testing  time. 


Table  10 

Successive  Approximations  to  I'D^  , for  = Tq  = 142 


Anprox*n 
(1)  I'D  : 


25.00  9.00  30.00  23.00  15.00  40.00  142.00 

25. 38  S.85  20.04  25.55  33.62  28.56  142.00 

25.83  7.60  16.02  25.41  42.43  24.73  142.02 

26.48  7.30  14.72  25.22  44.54  23.73  141.99 


Value  of  0 for 
Successive  Values 
of  L 

.265 
. 282 
. 283 


Table  11 

Successive  Approximations  to  I'D^  > Tor  = 2Tq  = 2(l42)  =■  284 


ADprox*  n 1 2 2 ± 5 6 Z 

(2)  I'D^i  1 50.00  18.00  60.00  46.00  30.00  80.00  284.00 

2 51.51  13.62  30.14  51.31  S3. 57  53-84  284.00 

3 54.84  10.84  22.51  51.17  97.72  46.91  284.00 

4 56.32  10.53  21.24  50.94  99.64  45.34  284.01 


Value  of  0 for 
Successive  Value 
of  L 

.305 

L2  .334 
I.J  .337 


t 
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III  MATHEMATICAL  DERIVATION 

In  (1)  a procedure  is  developed  fcr  altering  test  lengths  in  a battery 
to  give  maximum  multiple  correlation  vrith  a single  criterion.  The  develop- 
ment of  this  procedure  will  be  reviewed  and  the  procedure  will  be  extended  to 
the  problem  of  differential  prediction. 

Let 

M be  the  number  of  cases, 
n be  the  number  of  predictors, 

Z be  an  (M  X n)  matrix  of  test  scores  in  a battery  of  altered  lengths 

2 '*2 

with  the  elements  of  Z of  the  form  i.1  .1  , 

■/"m 

J 

W be  an  (M  x 1)  vector  of  criterion  scores  with  elements  of  the  form 
vr  - w 

v'li 

B be  an  (M  X 1)  vector  of  regre'sion  coefficients  for  estimating 
W from  Z, 

r be  an  (n  X n)  matrix  of  intercorrelations  of  tests  of  original  lengths, 

p be  an  (n  X n)  matrix  of  intercorrelations  of  tests  of  altered  lengths, 

r^  be  an  (n  X 1)  vector  of  validity  coefficients  for  the  tests  of 

original  lengths, 

be  an  (n  X 1)  vector  of  validity  coefficients  for  the  tests  of 
altered  lengths, 

D be  an  (n  X n)  diagonal  matrix  of  origins!  test  lengths, 

cl 

be  an  (n  X n)  diagonal  matrix  of  altered  tost  lengths. 

Dp  = be  the  ratio  of  altered  to  criginal  test  lengths, 

U be  the  (n  x n)  diagonal  matrix  of  reliability  coefficients  for  the 
^ii 

tests  of  original  lengths. 
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Let  6 = [l  - (V^)Dr..J^e  * 

Let  e = (ZB  - W)  . (2) 

We  wish  to  minimize  e'e  vrith  the  constraining  condi.tion  l^D^l  = '1 
where  T is  the  ratio  of  the  new  total  testing  time  to  the  original 
total  testing  time,  and  1 is  a column  vector  of  all  unit  elements. 

To  obtain  e'e  minimum  under  this  condition 
Let  = e'e  + M'D^I  (3) 

where  X is  a Lagrangion  multiplier. 

From  (2) 


^ = (B'z'ZB  - B'Z'W  - W'ZB  + W'W)  + Xl'D^l 
From  the  definitions  above 
Z'Z  = p 

z'w  = , 

W'V/  = 1 

Substituting  (5)>  (6),  and  (7}  m (4) 

it  = B'pB  - B'p  - p'B  ^ 1 + U'D,  1 . 

^ C D 

In  (l)  it  is  shown  that 


and 


P 


= 6”^(r  - D + D D 1X1)6“^ 

'■  u u a^' 

where  we  define  D “ I-D 
u r, 

unreliability  coefficients 
1 


a diagonal  matrix  of  test 


Let  B “ 6^3  . 

Substituting  (9)>  (10),  and  (11)  in  (8) 

= 3'  (r  - + D^D^D^1)3  - 3'r^  - r;3  + Xl'D^l  . 

The  unknowns  on  the  right  hand  side  of  (12)  are  3,  and  X . 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


(11) 


(12; 
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Differentiatinp,  (12)  v/ith  respect  to  p'  and  equating  tl:e  resulting 
expression  to  zero  to  get  an  extremum 

4'^  = r^-(r-D  + DD  D:^)p  = 0 

c ' u uab'*^ 

or  P = (r-D  + DD  (!' 

' u u a b'  c 

Differentiating  <jj  v/ith  respect  to  the  scalars,  b^,  (i  = 1,  2,  n) 

and  equating  the  n resulting  expressions  to  zero 

= 0 (li 


3i 


Sunuidng  these  n equations 

Thus  in  matrix  notation  v/e  obtain 

^ I'D  1 
b 

Substituting  for  in  (15)  and  collecting  these  n expressions  as  the 
diagonal  matrix  D^  we  obtain 

A I'V 

D,  = D„(D  D Y — I — 


where  D^  is  a diagonal  matrix  with  the  3^  as 
diagonal  elements. 


In  (1)  it  is  shown  that 


\ r - D + ^ 


I'd,  1 

D 
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Using  (18)  in  (1?)  we  can  therefore  solve  for  D^,  the  new  test  lengths. 

The  new  multiple  correlation  is  given  by 

- P'r,  . (19) 

Next  we  extend  the  procedure  to  the  case  of  differential  prediction. 
Consider  the  following  additions  to  the  definitions  given  above. 

Let 

N be  the  number  of  criteria, 

V/  be  an  (M  X N)  matrix  of  criterion  scores  whose  elements  are 
deviate  scores  of  the  form  i.1  .1  , 

3 

2 

H be  an  (M  X N ) matrix  consisting  of  difference  vectors  for  all 

possible  pairs  of  criterion  vectors  i and  j,  including  i =■  j, 

2 

B be  an  (n  X N ) matrix  of  "least  square"  regression  vectors  for 
estimating  H from  X, 

r be  the  (n  x N)  matrix  of  validity  coefficients  with  the  tests  of 
c 

original  lengths, 

p be  the  (n  x N)  matrix  of  validity  coefficients  vd.th  the  tests  of 
c 

altered  lengths. 

From  the  differential  prediction  procedure  (2)  we  have 

0 = I'D  1 - , the  index  of  differential  prediction  efficiency,  (20) 

C 


where  C = r'  r'^  r 
c c 

and  D is  a diagonal  matrix  of  the  diagonals  of  C, 
c 

(21) 

Let  E = ZB  - H 

(22) 

and  = e^l^  - I 

(23) 

where  e.  is  a column  vector  cf  all  zero  elements  except  the  i'th  which 
1 


is  unity. 
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Let 

G = (F^,  F2  F.^)  . 

(24) 

Thus 

we  have 

H = WG' 

(23) 

and 

E = ZB  - WG'  . 

(26) 

Fron  (23),  (24)  and  (26)  postmultiplied  by  C-  and  divided  by 

2N  we  obtain 

_ E G _ Z EG  W /l  - 11'  \ 

^ 2N  “ 2N  “ \ ‘Ij'*  I 

(27) 

since  G'G  = 2N  /l  - ll'  \ . 

1 N y 

Let 

(28) 

Let 

(29) 

Then 

e = ZJ  - t . 

(30) 

We  wish  to  minimize  the  trace  of  e^e  with  the  constraining  condition 

I'd,  1 * T . 

D 

Let 

= tr  e'e  + \1'D,1 

^ D 

(31) 

Let 

(32) 

Substituting  (5),  (6),  (9),  (10),  (28),  (29),  (30),  and  (32) 

in  (31) 

we  obtain 

- tr[j'6-'i 

U'D^l 

• (33) 

Let 

6"^J  = L,  > 

(34) 

Let 

1 I - u'' 

(35) 

Let 

li  = r - D 

u 

(36) 

Let 

A = D D 
u a 

(37) 

Substituting  (3L-),  (35),  (36)  and  (37)  in  (33)  we  obtain 

5^  = tr  ^_L'(R  + 4D^'^)L  - L'a^  - a'L  + t't  | + W.'Dj^1  . 

(3a) 
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Dif ferentiating  (3^)  with  respect  to  row  vectors  of  L'  and  equating  the 
results  to  0 w',  obtain 

fi,'  . 0 

or 

= (R  + AD^-i)L  . (39) 

Differentiating  (37)  with  respect  to  and  equating  the  results 
to  0 we  obtain 

^ - XI  - . 0 (40) 

u 

where  D, , < is  a diagonal  matrix  whose  non-zero 
elements  are  the  diagonal  elements  of  LL’  . 


(41) 

(42) 

(43) 

(44) 

(45) 


(46) 
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The  procedures  outlined  in  Section  II  may  be  related  to  the  above 
mathematical  development  as  follows; 

Table  1 is  given  by  equation  (36). 

Step  1 is  based  on  equation  (35), 

Step  2 is  based  on  equation  (37). 
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Step  3 is  based  on  equation  (46). 

Step  4 consists  of  oalculatiiig  from  . 

1 1 

Step  5 consists  of  calculating  . 

Step  6 consists  of  calculating  the  parenthesis  on  the  right  side  of 
equation  (45)  for  i = 1. 

Step  7 consists  of  calculating  from  equation  (45)* 

Step  8 consists  of  calculating  D.  from  equation  (47)« 

2 

Step  9 consists  of  calculating  nn  matrioc  from  equation  (45). 

Step  10  uses  equation  (47)  to  calculate  D.  . 

3 

In  general  steps  6 and  7 are  repeated  for  successive  vaiues  of  i in 
equation  (45)  and  step  8 is  repeated  for  successive  values  of  i in  equation 
(47). 

Step  11  uses  equation  (53)  to  get  successive  values  of  0. 
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